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Exotic stochastic processes are shown to emerge in the 
quantum evolution of complex systems. Using influence func- 
tion techniques, we consider the dynamics of a system coupled 
to a chaotic subsystem described through random matrix the- 
ory. We find that the reduced density matrix can display dy- 
namics given by Levy stable laws. The classical limit of these 
dynamics can be related to fractional kinetic equations. In 
particular we derive a fractional extension of Kramers equa- 
tion. 
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Whether one studies deterministic Hamiltonian or dis- 
sipative systems, one finds that transport in chaotic sys- 
tems often resembles some type of stochastic process. 
The dynamics of such systems leads to a rich spectrum of 
behaviors, from enhanced diffusion such as tracer diffu- 
sion in flows and turbulent diffusion in the atmosphere, 
to dispersive diffusion jlj]. Much effort has been spent 
in recent years to understand such classical stochastic 
processes in chaotic systems, leading to the development 
of approaches ranging from fractional kinetic equations 
Levy flights || to random walks in random envi- 
ronments HI! an d stochastic webs (?]] . One of the com- 
mon features to all of these is the use of Levy stable laws 
H . It was shown by Levy || , in studies of extensions of 
the central limit theorem, that a continuous class of non- 
gaussian processes satisfy the same fundamental equation 
that gives rise to the theory of gaussian processes, namely 
the Chapman-Kolmogorov equation for the conditional 
probability P(q, q', t): 



P(q,t/;t) = J dq"P(q,q",t-t")P(q 



" q'X). 



(1) 



solution, P(q,q',t) = cxp(— (q 



The standard 

q') 2 /4Z)i) / (4irDt) 3 / 2 , gives rise to the gaussian processes 
and the usual form of the Fokker-Planck equation. The 
general solutions of Levy provide a way to generalize 
Brownian motion. 

The non-gaussian processes which satisfy (1) are 
known as Levy stable laws, and have the form: 



P(q,t)=L£(q) 



1 

2^ 



exp{ifcg- A\k\ a }dk (2) 



to gaussian processes. The Levy distributions C^(q) sat- 
isfy the scaling relation: 



Ci(q) = A- 1 / a C 1 JqA- 1 ^), 



(3) 



where for A = 1 we drop the superscript: C\{x) = C a (x). 
For a < 2, these distributions are characterized by infi- 
nite second moments, as one can easily infer from the 
asymptotic behavior for q — > ±oo j^] , 
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These non-gaussian processes can be related to anoma- 
lous transport in a variety of (classical) physical systems 
||, as well as to classically chaotic systems. We have re- 
cently shown that turbulent diffusion can also arise in the 
time evolution of complex quantum systems [ [Tof . Here 
we find that a general form of quantum chaotic back- 
grounds can give rise to quantum evolution characterized 
by Levy distributions. Further, we can connect, in the 
semi-classical limit, such processes to fractional kinetic 
theory, which was initially introduced as a phenomcno- 
logical approach to classical anomalous diffusion. 

We would like to study the problem of a particle cou- 
pled to a chaotic environment, quantum mechanically. 
It has been realized in recent years that the quantum 
counterpart of chaos is random matrix theory. For sys- 
tems with time-reversal symmetry, the random matrices 
are real-symmetric. In this Letter we will examine the 
class of quantum dynamic processes, which can be real- 
ized through the interaction of a particle with a random 
matrix background. In contrast to the Caldeira-Leggett 
approach , we assume from the outset that the back- 
ground is chaotic, and not necessarily thermal. We de- 
note the coordinates of the background by (x,p) and that 
of the test particle by (X, P) . The Hamiltonian for the 
background plus interaction is taken to have the following 
form: 

H b = h (x,p) + hx(X,x,p). (5) 

In the basis of (many-body) eigenstates of ho, ho \ n) — 
e n | n) (n = 1, N), we define the matrix of Hb as 



[H b ]ij =e l 8 lj + [hi(X)] 



(6) 



where < a < 2 and A cx t. The case a = 2 corresponds 



It is convenient for calculations to choose an average level 
density as p(e) = poexp(/3e). For a background with 



1 



constant average level density, (3 = 0, while for a general 
many body system, j3 > 0. The chaotic properties of the 
interaction of the background with the particle are incor- 
porated into the correlation function (second cumulant): 



(([/n P0k#i OHM = Ga{x - Y)A ijkl . 



(J) 



Here Ajjjy = [5^6 ji + Sudjk], and all other cumulants 
vanish. In our analysis, the integration over the chaotic 
part, given by h\(X), is defined through a gaussian mea- 
sure for parametric random matrices pig] 

P[hi(X)] oc exp j-~ y dAdY~ (8) 

Tr[/ ll (X)a- 1 (X-y))/i 1 (Y)]}. 

The character of the interaction of the background with 
the test particle is incorporated into the correlation func- 
tion Q(X — Y), for which we use the form 



Sij(X) 



2Try/p(e l )p(e 



■. exp 



2kq 



G 



(9) 



This describes a parametric, banded, random matrix 
where the strength of matrix elements decreases with in- 
creasing level density. Here G(x) = G(—x) — G*{x) < 1, 
G(0) = 1, and the spreading width , n (linked with 
the effective band width N K p(e)) and the correla- 
tion length Xq are characteristics of the background. 

In order for the measure (||) to be positive definite G 
must not decorrelate faster than a gaussian fl2]| : 



G(x) 



a e (0,2]. 



(10) 



As the position X of the slow particle changes, the instan- 
taneous energy levels E n (X) of [Hb(X)]ij change. Using 
the above measure, the average fluctuations are 



([E n (X) - E n (Y)} 2 ) =D a \X- Y\ a . 



(11) 



The energy-spacing fluctuations have a behavior, which 
is similar to a Levy process characterized by the diffusion 
constant D a . The character of these fluctuations in the 
eigenvalues E n (X), indicated by a, will be seen to be 
related to Levy distributions, which describe the time 
evolution of the density matrix for a particle evolving in 
this chaotic bath. 

To develop the dynamical evolution of a free particle 
evolving in the presence of a chaotic background, we take 
the Hamiltonian of the form: 



Hij(X,P) 



p2 

2M+™ 



+ H b<ij {X). (12) 



The correlated, random-matrix bath can be integrated 
out in an influence functional formalism J3. For our 
purposes, the o(/3) action is sufficient, as well as a weak 



coupling of the particle to the bath. In this case the ef- 
fective equation for the density matrix of the test particle 
has the form: 

ihd t p(X, Y, t) = { H - |i + U(X) U(Y) (13) 



4XnM 
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X< 



(Px - Py) 
p(X,Y,t), 



where in weak coupling, G(x) = 1— |a;| Q and G'(x) above 
represents — asign(x)|a;| Q_1 . 

Consider first a test-particle in the absence of an ex- 
ternal field and interacting with a background with con- 
stant average level density (U(X) — and (3 = 0). This 
evolution equation can be solved by passing to the coor- 
dinates r = (X + X')/2, s = X — X' . In these variables, 
the density matrix has the form: 



p(r,s,t) 



dr 

T l M 
Tik 



dk 



kt 



— po \ r , s exp 

2nn \ m 



f ds' 












Js-kt/M 







ik(r — r') 



(14) 



An initial wavepacket, iJ}q(X) = 

cxp[— A 2 /4cr 2 ]/[27rcr 2 ] 1 / 4 , provides an initial density ma- 
trix Po (X,X') = (l/V^W)expH4r 2 + s 2 )/8a 2 ]. 

For the diffusive dynamics of the test particle, we are 
interested in the diagonal component of the density ma- 
trix p(X, X, t) = p(r, s = 0, t) : 



p(r,0,t) 



dr'dk 
2nh 



po (r',~fci/M)exp 



,MT l 



ik(r — r') 



ds 1 

kt/M Kn 

dk 
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(a + l)h(MX ) a ^ a+lk h 



(15) 



(16) 



p(X, X, t) is nothing more than the spatial probability 
distribution P(X, t) for the process. We can now express 
it in terms of a convolution of Levy distributions: 



p(X, X, t) = J dX'C a a {t) {X')C b 2 (t) (X - X') 



where 



a(t) = 



ri ( n 



(a + l)h \MX C 



ia+1 



b(t) 
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(17) 

(18) 
(19) 
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As both functions in the integrand of Eq. (17) are pos- 
itive definite, the spatial probability P(X,t) is also pos- 
itive definite. Notice that the restriction of < a < 2, 
which came from the short-distance statistical correla- 
tions in (11) and the requirement of a positive definite 
statistical measure, is also the necessary requirement on 
the Levy distribution to keep the resulting time evolu- 
tion positive definite. Hence the character of the short- 
distance fluctuations is directly responsible for the long- 
time behavior of the quantum system. 

Consider now the short-time and long-time behavior 
of the dynamics. For 1 < a < 2, in the limit of long times, 
we expect the t a+1 term to dominate over t 2 in (|l6|), so 
that the density asymptotically approaches a Levy dis- 
tribution: 



p(X,X,t)— » a(t)-^ a C a (a(i)- 1/Q x) , 



(20) 



while for very short times, the gaussian process is the 
dominant behavior: 



p(X,X,t) 



— £2 { — X 



(21) 



Specifically, in Eq. ( jig ) the \k\ a term in the exponent 
dominates in the long time limit only for momenta k < 
k c , where 



k r . 



b(t) 



oc t 2 - 



(22) 



For the special case of a = 2, the result is gaussian, 
but the dynamics can be anomalous, as one can have 
turbulent-like diffusion of the type (X 2 ) ~ t 3 @. When 
the level density of the background is not constant, (3 > 
and a = 2, one can recover Brownian diffusion Jl5| . 
For general a and f3 > 0, however, the results are not 
yet known. For the range < a < 1, the long-time 
behavior approaches a gaussian process. At short-times, 
the dynamics is influenced by C a , and there is a cross- 
over from short time Levy dynamics to normal gaussian 
expansion of the wavepacket. In both cases however, 1 < 
a < 2 and < a < 1 respectively, the second spatial 
moments are strictly speaking divergent. One should also 
note that as neither a(t) nor b(t) are linear in time, even 
though the dynamics has the character of a Levy process, 
it is not a Levy stable law. 

Efforts to understand unusual stochastic behaviors of 
dynamical systems has led to the development of exten- 
sions of the Fokker-Planck (FP) equation |-§. These 
are phenomenological fractional kinetic equations (re- 
stricted to one dimension) in which certain derivatives 
are replaced by derivatives of 'fractional' order . Such 
approaches have also found applications in a wide range 
of problems from turbulence to diffusion in porous or 
viscoelastic media [pi . We can now explore the type of 



stochastic process, which emerges in the classical limit 
of our quantum Levy processes, and the connection to 
multi-dimensional fractional kinetic theory. 

A typical type of phenomenological fractional FP equa- 
tion has the form 



d s P(Q,t) 
dtP 



(A(Q)P(Q,t)) 



d{-QY 



(23) 



where p = v — 1 in Ref. Q, p = v in Ref. ||] and 
8 = v = 1 in Ref. Q. Here the symbol d^/dx^ 1 rep- 
resents the Riemann-Liouville fractional derivative (jl7j, 
except for Ref. || , where it represents the Fourier trans- 
form of — fc M . This equation, while formally constructed, 
is phenomenological. It is defined to reproduce anoma- 
lous diffusion through scaling formulas such as Q 2 ~ t 1 , 
where 7 is a function of <5, p, v. A few points should 
be made here. Generally, the coefficients A and B are 
defined as limits whose existence is postulated but not 
known. Further, either the form of the fractional deriva- 
tives is taken to provide this scaling law, or power law 
noise is chosen to obtain them ^Jj|. Such dynamics can 
then be related to Levy processes ||]. Finally, the ex- 
tension of these equations to phase space becomes tenu- 
ous, since it is not clear how to include momentum. Not 
only is it unclear if one should take fractional derivatives 
with respect to coordinates, momenta or both, but the 
existence of the corresponding coefficients A, B,... is un- 
known. Through our transport equation, we can provide 
a microscopic interpretation of these coefficients as well 
as a systematic manner to construct a fractional kinetic 
equation in phase space, whose quantum limit results in 
Levy processes. 

To obtain a classical transport equation, we construct 
the Wigner transform f(Q,P,t) of the density matrix 
p(X, Y, t) as 



/ = 



dR 
2nh 



exp 



R 



,Q 



R 



,t . (24) 



Applying this to our evolution equation, taking the lead- 
ing order terms in H, we find 



df 

at 



dR 



■ exp 



iPR 



2M 



d Q d R 



(25) 



R 



X 



+i r jhXQasigii(R) 



R 



Xn 



9r > P 



R „ R 



This leads naturally to the Reisz fractional integro- 
differential operator. This operator, applied to a function 
f(P), is defined as @ 



(-A P )*/ = T- 1 \X\ a Tf 1 



(26) 
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where Ap is the Laplacian (in our case with respect to 
the momentum P) , and T represents a Fourier transform. 
(This operator is distinct from that proposed in which 
did not have the absolute value, and from ||, which uses 
the Riemann-Liouville form of this operator. The Reisz 
operator is defined as a fractional integral for Re a <0 
and as a fractional derivative for Re a >0 through ana- 
lytic continuation.) It is convenient to define the operator 
Dp = {~i/h) a (~A P ) a / 2 , since D%[f] = d 2 f/dP 2 . Then 
the classical limit of our quantum Levy process gives rise 
to a fractional extension of Kramers equation: 

df(Q,P,t) P df(Q,P,t) dU(Q)df(Q,P,t) 

dt M dQ dQ OP { ' 



Levy dynamics. In passing to the classical limit of this 
behavior, we find that the dynamical evolution results 
in a fractional kinetic equation, which is a generaliza- 
tion of Kramers equation. For a — 2 Kramers theory is 
recovered. This approach provides a means to develop 
fractional kinetic theory in more than one dimension, 
since the expansion coefficients are determined from the 
microscopic theory. It also provides the possibility to 
explore the connections between quantum and classical 
transport in chaotic systems, as well as the links between 
chaos, quantum statistical fluctuations, Levy processes 
and classical fractional dynamics. 



la {D a p [Pf(Q,P,t)} 



2TM 
ah 2 



(iK) a D a P [f{Q,P,t)} 



where T = 1/ ft is the temperature and the operator 
D P = (-i/h) a J r - 1 sign(X)\X\ a F, with the property 

D P [Pf] = d(Pf)/dP. The generalized friction coeffi- 
cient is given by: 



7o 



2M Xn ' 



(28) 



For a — 2 we recover Kramers equation jL8||. What we 
see is that it is not the coordinates which acquire the frac- 
tional character, as usually assumed, but the momenta. 
Because the coupling to the background is not momen- 
tum dependent, the correlation function G(X) results 
only in fractional derivatives with respect to momenta. 
This can be traced back to the nature of the chaotic 
correlations in Eq. (13). Further, these processes, re- 



lated to Levy processes, do not require the introduction 
of fractional time derivatives. We note here that our 
transport theory has a consistent classical limit for all of 
these transport coefficients only when they remain finite 
as h — > 0. This requires in turn that the parameters of 
our quantum theory cannot remain constant as h — > 0, 
if we are to recover a well defined classical transport. 
Finally, we observe that this approach provides finite co- 
efficients Dqq, Dpp, Dqp and so forth (e.g. A, B,...) 
for a fractional kinetic equation in phase space. 

We have shown that the quantum evolution of a 
wavepacket in a chaotic environment can lead to reduced 
density matrices which behave as Levy processes. The 
short distance energy fluctuations of the background, 
which are characterized by a parameter a G (0,2], are 
found to be precisely related to the quantum time evo- 
lution with a Levy process of the same character a. 
For a = 2 one has gaussian processes, which can dis- 
play normal to turbulent -like diffusion or Brownian dif- 
fusion (/3 > 0), while for a — 1 one has the dynamics of 
the Dyson process. The general quantum evolution of a 
wavepacket displays a cross-over between gaussian and 
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